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Abstract 

We report the results of a search for the rare baryonic decay modes B° — ► 
AA, B + — > Ap, B° — > Ap7r~, and 5° — ► (and their charge conjugate 
states) using 5.8 x 10 6 BB pairs collected with the CLEO detector. We see no 
statistically significant signals in any of these modes and set 90% confidence 
level upper limits on their branching fractions, B(B° — > AA) < 3.9 x 10 -6 , 
B(B+ Ap) < 2.6 x 10~ 6 , B(B° -» ~Kjm~) < 1.3 x 10~ 5 , and 5(5° -> pp) 
< 7.0 x 10~ 6 . 
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FIG. 1. Tree (a,c,e) and penguin (b,d,f) decay processes which are expected to dominate the 
baryonic decays under consideration. 



Evidence for the b — > s quark transition, allowed in the Standard Model by the pen- 
guin (internal W boson-quark loop) diagram, has been observed in both electromagnetic 
decays and hadronic decays of B mesons to two mesons ||f§]. Charmless B meson de- 
cays can also arise through the b — > u tree process as well as the Cabbibo suppressed b — > d 
penguin process. These rare decay processes are of considerable theoretical and experimen- 
tal interest due to their importance in probing the Cabibbo-Kobayashi-Maskawa (CKM) || 
picture of CP violation within the Standard Model. 

This communication describes a search for charmless baryonic B decays to the final 
states AA, Ap, Apir~ , and pp (and their charge conjugate states). The dominant tree-level 
and one-loop penguin diagrams contributing to these decays are shown in Figure |l| The weak 
processes are similar to the meson decays B — > tttt, Ktc, and KK. We expect B + — > Ap and 
B° — > Apn~ to be dominated by the b —>■ s penguin diagram, while B° — > AA and B° — ► pp 
should be dominated by the b — ► u tree process. W-exchange, annihilation, penguin annihi- 
lation, and electroweak penguin processes can also contribute to these decay amplitudes ||. 
However, these diagrams are expected to have small contributions as compared to those 
shown in Figure [TJ. Theoretical predictions for these modes, scaled for a common value of 
\Vub\ = 0.0033 J7j, are given in Table |T[ Both pole model || and QCD sum rule calcula- 
tions |J have been performed. The pole model calculations yield higher predicted branching 
fractions for the rare baryonic modes. The search for the three-body final state Apn~ is 
motivated by the experimental observation that B(B~ — > A+p7r~)/£>(£>° — > A+p) > 1.0 
(90% CL) [|H]]. We estimate the branching fraction for Apir~ from the measured branching 
fraction of B~ — > A+pvr - , scaling by the appropriate CKM and phase space factors. Penguin 
contributions could lead to significant enhancement over this rough estimate. 



The data sets used in this analysis were collected with the CLEO II [11] detector at the 
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Cornell Electron Storage Ring (CESR). It consists of 5.41 fb~ taken on the T(4S) resonance 
and 2.79 fb _1 taken below the BB production threshold. The on- resonance sample contains 
3.3 million BB pairs taken before the installation of the silicon vertex detector, |L2] (CLEO 
II data set) with the remainder of the 5.8 million BB pairs taken after the detector upgrade 
(CLEO II. V data set). 

The momenta of charged particles are measured in a 67-layer tracking system operating 
inside a 1.5 T superconducting solenoid. The main drift chamber also provides a measure- 
ment of the specific ionization loss, dE/dx, which is used for particle identification. Photons 
are detected by the 7800-crystal Csl calorimeter. Muons are identified using proportional 
counters placed at various depths in the steel return yoke of the magnet. 

Charged tracks are required to pass track quality requirements based on the average hit 
residual and the impact parameters in both the r — <fi and r — z planes. Pairs of tracks 
with vertices displaced by at least 5 mm from the primary interaction point are taken as 
A candidates. We require the pir~ invariant mass to be within 10 MeV/c 2 of the A mass. 
Reconstructed A baryons have a mass resolution of 1.4 MeV/c 2 . The impact parameter and 
A flight distance requirements nearly eliminate feed-across between the AA and Ap-K~ final 
states. 

Charged particles are identified as protons or pions from specific ionization (dE/dx) 
measurements from the drift chamber. Electrons are rejected based on their dE/dx and the 
ratio of the track momentum to the associated shower energy in the Csl calorimeter. We 
reject muons by requiring that the tracks do not penetrate the steel absorber to a depth 
greater than five nuclear interaction lengths. For the pp final state, a requirement on the 
time-of-flight of the protons and antiprotons is applied that eliminates 40% of pion and 28% 
of kaon background while retaining 89% of the signal. 

We calculate a beam-constrained B mass Mbc = yE~l ~~ Ps^ where p_e is the B candidate 
momentum and Ey, is the beam energy. The Mbc resolution is about 2.5 MeV/c 2 for all 
modes. We can capitalize on energy conservation if we define AE = J2i{Ei) — E\,, where Ei 
are the energies of the daughters that form the B meson candidate. The AE distribution is 
then centered about zero for reconstructed B mesons and has a mode-dependent Gaussian 
width that ranges from 17 MeV (for Apix~) to 25 MeV (for pp). We accept events with 
M BC within 5.2 — 5.3 GeV/c 2 and \ AE\ < 200 MeV. This region includes signal as well as 
a sideband used to fix the background normalization. 

Backgrounds from b — > c as well as other b — > u and b — > s decays are negligible for 
the two-body decay modes, since the signal daughter particles are relatively light and are 
produced with high momentum. This is in contrast to the dominant B meson decays which 
typically have a large number of final state particles with lower momenta. The decays 
B — > K + tt~ and B° — > tt + tt~ are displaced from zero by +287 and +330 MeV in AE, 
respectively, and do not pass the sideband AE requirements for the pp final state. Possible 
backgrounds from B decays to pseudoscalar-vector final states (B — > PV) such as B — > p-n 
and B — > K*ir are found in simulation studies to contribute at most 4 ± 4 events to pp. The 
simulations are normalized to 90% CL upper limits for B{B — > PV) |R|. For the three-body 
mode Apir~ , we suppress backgrounds from B meson decays as well as continuum by requiring 
that the p, ir~ , and A momentum be greater than 0.70, 0.75, and 1.0 GeV/c, respectively. The 
placement of these cuts was determined by a signal squared over background optimization. 
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TABLE I. Experimental results and theoretical predictions [8,9]. Branching fraction (B) upper 
limits at 90% CL are given in 10 -6 units with (without) systematic errors included. Quoted 



significance of the fit result is statistical only. 

Mode s Sig. £(%) B Theory B 

AA O-Olao °- 0cr 14.8 ±1.8 < 3.9 (3.4) 0.13 

pp mtti 2 - 8cr 48.8 ±2.7 < 7.0 (6.0) 0.6-4.8 

Ap O.Ol£o O.Ocr 29.3 ± 2.8 < 2.6 (2.3) ~ 3 

Apn- 3.0±|;| l.ler 14.6 ±1.7 < 13 (11) -0.5 



These requirements reduce the background from other B decays to less than 5% of the total 
background. 

For this analysis, the main background arises from e + e" — > qq (where q = u,d,s,c). 
Such events typically exhibit a two-jet structure and produce high momentum back-to-back 
tracks. To reduce contamination from these events, we calculate the angle #t between the 
thrust axis of the candidate tracks and the thrust axis of the tracks and showers in the rest 
of the event. The distribution of cos#t is strongly peaked at ±1 for qq events and is nearly 
flat for BB events. For the low background final state AA, we require | cos(9t| < 0.95 and 
for the remaining modes we require | cos^t] < 0.9. A cut of | cos#t| < 0.9 eliminates 66% 
of continuum background while retaining 87% of the signal. 

A detailed GEANT-based Monte Carlo (MC) simulation |14 was used to determine the 
overall detection efficiencies {£) for each mode as shown in Table |I[ Efficiencies contain 
branching fractions for A — > where applicable. For the efficiency calculation we assume 
that the decay B° — > Ap-K~ proceeds by phase space. |1| From independent Monte Carlo 
samples we estimate systematic errors on the efficiency determination; those errors are also 
given in Table |I[ 

Additional discrimination between signal and qq background is provided by a Fisher 
discriminant technique, described in detail in Ref. [|Lj. The Fisher discriminant, T = 
J2}Li a iVi-, is a linear combination of 11 variables (yi) where the coefficients (c^) are chosen 
to maximize the separation between signal and background Monte Carlo samples. The input 
variables are | cos 6^1 (the cosine of the angle between the candidate thrust axis and the beam 
axis), | cos#b| (the cosine of the angle between the B meson momentum and the beam axis), 
and the energy and momentum contained within nine concentric 10° cones that surround 
the candidate thrust axis. The sum of energy and momenta from tracks and showers in the 
forward and backward cones are combined. The statistical separation between signal and 
continuum background afforded by the Fisher discriminant is 1.5 standard deviations (a) 
after events with | cos^tI > 0.9 have been rejected. 

To determine the signal yield in each mode we perform an unbinned maximum likelihood 
(ML) fit using AE, Mbc, A and A masses (where applicable), and dE/dx as input 
information for each candidate event. Separate fits are performed for each mode. The 
likelihood function is defined as 
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N 



£ = J] [sPsixi) + bP B (%)] 



i=l 



where s (b) are the signal (background) yields in the candidate sample (of N total events), 
Ps and Pb are the signal and background probability density functions (PDF), respectively, 
and Xt are the appropriate input variables discussed above. The signal and background yields 
are varied until the likelihood function is maximized. The probability density functions Ps 
and Pb are formed by the product of the best fit functional forms for signal and background 
distributions for each input variable. Signal MC is used to determine the shape of the signal 
events and off-resonance data is used to determine the background shapes. We use the sum 
of two Gaussians to parameterize the signal shapes for AE, Mbc, an d A masses. Bifurcated 
Gaussians (with different low and high side widths) best describe the signal and background 
shapes for the Fisher discriminant and a straight line is used to parameterize the AE back- 
ground shape and non-resonant A mass background. When determining the background 
Mbc distribution |16] the kinematic endpoint of the off-resonance data is matched to the 
on-resonance data by shifting the mass distribution by 5.290 GeV — E^. 

Systematic errors on the fitted yields are determined by varying the PDFs used in the fit 
by their measurement errors. This variation is performed with a Monte Carlo technique in 
which the PDF parameters are varied by Gaussian distributed random numbers. The width 
of the Gaussian is equal to the measurement error of each parameter. For each variation 
the data sample is refit. The distribution of the fit yield and upper limit yield from these 
MC experiments is fit to a Gaussian whose width determines the ML fit systematic error. 
Important correlations between PDF parameters are accounted for in this technique. 

For each of the four rare baryonic decay modes, Table | shows the signal yield (s) and 
statistical significance of the yield as determined from the likelihood function. Since we see 
no statistically significant signals in any of the modes, we calculate 90% confidence level 
upper limit yields (s UL ) by integrating the maximized likelihood function (£ max ) as follows: 

f °° L max (N)dN " U - yU - {Z) 

The upper limit yield is then increased by its systematic error (from the PDF variation) and 
the detection efficiency reduced by its systematic error before we calculate the branching 
fraction upper limits shown in Table |I[ Table [I] also gives the upper limits for each mode 
before systematic errors have been included. 

Figures |2] and |3| show the M B c and AE projections, respectively, for each of the four 
modes. To reduce the background level in the plots we require (cos^t) < 0.8 and apply 
a Fisher discriminant requirement, which eliminates roughly 80% of the background and 
retains about 60% of the signal. For the higher background modes (Ap, Apn~, and pp) we 
require |Ai?| < 2.5 (Tae for events to enter the M B c plot and \M B c — 5.28| < 2.5 ct Mbc 
for events to enter the AE plot. Signal efficiency with these stricter requirements is about 
50% of that quoted in Table p] while background efficiency is only 1-3 %. Overlaid on these 
plots are the projections of the PDFs used in the fit, normalized according to the fit results 
multiplied by the efficiency of the additional constraints. Both signal and signal upper limit 
yields are illustrated in the figures. 
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FIG. 2. M BC projection plots for (a) B° -► AA, (b) 5+ -► Ap, (c) B° -► Ap7r~, and (d) 
5° — > pp. The scaled projection of the total likelihood fit (solid curve) and the fit reported upper 
limit (dotted curve) are overlaid. 
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FIG. 3. AE projection plots for (a) B° AA, (b) 5+ -► Ap, (c) B° -► Apvr", and (d) B° 



pp. 



The scaled projection of the total likelihood fit (solid curve) and the fit reported upper limit (dotted 
curve) are overlaid. 
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We have searched for decays of B mesons to the baryonic final states AA, Ap, Ap7t~, and 
pp (and charge conjugate states) in 5.8 million BB pairs collected with the CLEO detector. 
We see no statistically significant evidence for signals in any of these modes and set upper 
limits on their branching fractions. The upper limit for B + — > Ap is slightly lower than the 
theoretical estimate given in Reference |J . We see no evidence for b —>■ s penguin transitions 
(recently observed in decays to meson final states) in B decays to baryonic final states. 
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von Dombrowski thanks the Swiss National Science Foundation, and H. Schwarthoff thanks 
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tional Science Foundation, the U.S. Department of Energy, and the Natural Sciences and 
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